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Abstract . In an expanding universe the vacuum energy density p\ is expected to 
<]} be a dynamical quantity. In quantum field theory in curved space-time p\ should 
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exhibit a slow evolution, determined by the expansion rate of the universe H. Recent 
measurements on the time variation of the fine structure constant and of the proton- 
electron mass ratio suggest that basic quantities of the Standard Model, such as the 
QCD scale parameter Aqcdj may not be conserved in the course of the cosmological 
evolution. The masses of the nucleons tun and of the atomic nuclei would also be 
affected. Matter is not conserved in such a universe. These measurements can be 
interpreted as a leakage of matter into vacuum or vice versa. We point out that the 
amount of leakage necessary to explain the measured value of m,N /mjsr could be of the 
same order of magnitude as the observationally allowed value of pa/ PA, with a possible 
contribution from the dark matter particles. The dark energy in our universe could be 
the dynamical vacuum energy in interaction with ordinary baryonic matter as well as 
with dark matter. 



PACS numbers: 95.36.+X, 04.62. +v, ll.10.Hi 
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1 Introduction 



The Standard Model (SM) of the strong and electroweak (EW) interactions contains 27 indepen- 
dent fundamental constants: the QED fine structure constant a cm = e 2 /47r, the SU(2)l gauge 
coupling g of the EW interactions, the gauge coupling constant of the strong interactions g s , the 
mass M\y of the weak gauge boson W, the mass M% of the Higgs boson T~L, the 12 masses of the 
quarks and leptons, the 3 mixing angles of the quark mass matrix, a CP-violating phase, the 3 
mixing angles in the lepton sector, a CP- violating phase and two additional phases, if the neutrino 
masses are Majorana masses. One of the parameters in the list, the mass of the Higgs boson M%, 
has not been measured thus far, despite some recent hints pQ. 

If we include the Einstein-Hilbert (EH) Lagrangian of gravity, there are two more fundamental 
constants, both of them dimensionful: Newton's gravitational coupling Gn and the cosmological 
constant A (also denoted as the CC term). The gravity constant has the dimension of an inverse 
mass squared (in natural units): G = 1/Mp, where Mp ~ 1.22 x 10 19 GeV is the Planck mass, the 
largest mass scale in the universe. The cosmological constant has the dimension of mass squared, 
the mass being of order Hq ~ 10 -33 eV, i.e. essentially the value of the Hubble parameter at 
present (the smallest mass scale in the universe). 

Until recently the observational data on |A| could only place an upper bound, but now cosmo- 
logical observations give a value, which is tiny, but non-vanishing (in particle physics standards) 
and positive [21 [3] . It can be expressed as an energy density: p\ = A/(8ir Gn) ~ 10~ 47 GeV 4 - 
the so-called vacuum energy density. We can define the mass scale associated to the CC term 
as follows: ttia = (p^) 1 ^ ~ 10 -3 eV. The scale m\ is the geometric mean of the two extreme 

1 /2 

mass scales in the universe: m\ ~ (Hq Mp) ' . In the ACDM model (i.e. the standard model of 
cosmology) this scale associated to the vacuum is assumed to be constant. This is a big puzzle 
within the ACDM model. 

The dark energy (DE) problem was originally presented as the cosmological constant (CC) 
problem [H El E] • This is one the basic problems of physics, ever since it was first formulated 45 
years ago [7]q 

In this paper we suggest the possibility that some of the cosmological constant problems might 
be related to basic parameters of the Standard Model. The nucleon mass and the QCD scale 
Aqcd might not have remained constant throughout the history of the cosmological evolution [9l 
[TUl \TT\ 112] , This is related to the time variation of the fine structure constant. Constraints on 
the ratio d cm /a cm can be derived from limits on the position of nuclear resonances in natural 
fission reactors which have been working for the last few billions years - the so-called "Oklo 
phenomenon" [131 O IS] ■ 

There could also be a cosmic time variation of the strong coupling constant, a s , related to the 
variation of the fundamental QCD scale Aqcd- One expects that Aqcd /Aqcd should be larger 
than a cm /a em . Recent high precision experiments performed both in the laboratory with atomic 
clocks \T7\ [T8] and in astrophysics using data from quasars [19] support these ideas. 

It has been suggested that the parameters of the EH action, Gn and A, may be varying with 

x For a recent detailed account of the old fine tuning CC problem, see e.g. sect. 2 and Appendix B of [jjj- 
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time due to the interaction of the vacuum with the mattero This time evolution might be linked 
to the time variation of the QCD scale and to the time shift of all the particle masses, including 
the dark matter ones. All the atomic masses of the chemical elements would be affected. 

Here is the outline of this paper. In section [2] we review the models with time evolving cosmo- 
logical parameters. In section[3]we specialize to a class of these models, where the time evolution is 
viewed as a renormalization group evolution. In section|4]we describe some experiments, providing 
evidence of the time variation of masses and couplings, and suggests a link of this variation with 
that of the cosmological parameters. In section [5] we propose the existence of a leakage of matter 
into the vacuum as a possible source of dynamical dark energy and compute the variation of the 
particle masses and the QCD scale with the Hubble rate. The last section contains our conclusions. 

2 Cosmological models with time evolving parameters 

We discuss the possibility that the cosmic time variations of the constants of particle physics and 
of cosmology are related. Scenarios, in which G could be variable, have been previously discussed 
in the literature. Dirac suggested in the thirties (through his "large number hypothesis" [21])) that 
the gravitational constant G could be varying with time in correlation with other fundamental 
constants. We also mention the ideas on time varying fundamental constants by Milne and Jordan 
at about the same time |22] . Later the time variation of G was tied to the existence of a dynamical 
scalar field coupled to the curvature - the original Jordan and Brans-Dicke proposals [23] . 
Consider the General Relativity field equations in the presence of the cosmological term: 

G^ u - g flu A = 8nGT flu . (2.1) 

Here G^ u = R^ u — ^g^R is the Einstein tensor, and is the energy-momentum tensor of the 
isotropic matter and radiation in the universe. Without violating the Cosmological Principle within 
the context of the FLRW (Friedmann-Lemaitre-Robertson- Walker) cosmology, nothing prevents 
the parameters G = G(t) and A = A(t) to be functions of the cosmic time, as it is the case with 
the scale factor itself a = a(t). The possibility of a variable CC term has been considered by many 
authors from different points of view [244 ES] , including the more recent quintessence approach - 
cf. [5] and references therein. 

The contribution from the A term, originally on the l.h.s. of Einstein's equations, can be 
absorbed on the r.h.s. after introducing the quantity p\ = A/(8"7rG/v), which represents the 
vacuum energy density associated to the cosmological term. Einstein's equations can then be 
rewritten formally the same way as in (|2.ip . but replacing the ordinary energy- momentum tensor 
of matter by the total energy-momentum tensor of matter and the vacuum energy: 

Tfj, u -> T^ v = T^y + g^ v p A = (pa ~ Pm) 9p,v + (Pm +Pm)U fJ/ U u . (2.2) 

Here p m andp m are the proper density and pressure of the isotropic matter, and is the 4- velocity 
of the cosmic fluid. 

2 For a recent review, see e.g. [20] and references therein. 
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The corresponding equation of state (EoS) oj m = Pm/Pm reads: uj m = 1/3 and uj m = 0, for 
relativistic and non-relativistic matter respectively. The redefinition of the energy-momentum 
tensor can be done in the same way as in Eq. (|2.2p . whether p\ is strictly constant or time varying. 
In both cases it enters with the equation of state p\ = —p\, i-e. u}\ = — 1. This is in distinction 
to the general DE fluids, whose EoS take the generic form pjj = ojdPd (with wo < —1/3) [21 E]. 

We discuss now some possible scenarios for variable cosmological parameters that appear when 
we solve Einstein's equations (|2.ip in the spatially flat FLRW metric, ds 2 = dt 2 — a 2 (t)(ix 2 , where 
a(t) is the time-evolving scale factor. We restrict ourselves to the spatially flat case, since this 
seems to be the most plausible possibility in view of the present observational data [2] and the 
natural expectation from the inflationary universe. We consider Friedmann's equation with non- 
vanishing p\, which provides Hubble's expansion rate H = a/a (d = da/dt) as a function of the 
matter and vacuum energy densities: 

H 2 = ^-(p m+ p A ). (2.3) 

As stated, we assume that pa = PA(t) and G = G(t) can be functions of the cosmic time t. We 
will denote the current value of the Hubble rate by Ho = 100 hKm/s/Mpc. The observations give 
h c± 0.70. The dynamical equation for the acceleration of the universe is: 

a 47rG, n n . 4ttG „ , 8tt G , n AS 

- = — {Pm. + 3p m - 2p A ) = — (1 + 3u m ) p m H — /9A • (2.4) 

a 3 3 3 

In the late universe (p m — > 0) the vacuum energy density pa dominates. It accelerates the cosmos 
for pa > 0. This may occur either, because pa is constant, and for a sufficiently old universe one 
finally has p m (t) < 2pA, or because pa(£) evolves with time, and the situation p\(t) > p m (t)/2 
is eventually reached sooner or later than expected. The general Bianchi identity S/^G^ = 0, 
involving the Einstein tensor on the l.h.s. of Eq. f|2. 1 [) . leads to the following relation for the full 
source tensor on its r.h.s. (after we include the CC term): 

V" (G T^) = V M [G (T^ + 9llv p A )] = . (2.5) 
The last equation provides the following "mixed" local conservation law: 

j t [G(p m + p A )]+3GH(p m +p m ) = 0, (2.6) 

where G and/or pa may be functions of the cosmic time. Although the previous equation is not 
independent of (|2.3p and (|2.4p . it is useful to understand the possible transfer of energy between the 
vacuum and matter, with or without the participation of a time-evolving gravitational coupling. 
For instance, if pA ^ 0, matter is not generally conserved, since the vacuum could decay into 
matter, or matter could disappear into vacuum energy (including a possible contribution from a 
variable G, if G ^ 0). The local conservation law (|2.6p mixes the matter-radiation energy density 
with the vacuum energy pa- 

We mention the following possibilities: 

• Model I: G =const. and pa =const.: 
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If there are no other components in the cosmic fluid, this is the standard case of ACDM 
cosmology, implying the local covariant conservation law of matter-radiation: 



p m + 3 H (p m + Pm ) = 0. (2.7) 



• Model II: G =const and p\ ^ 0: 

Here Eq. ([2.6p leads to the mixed conservation law: 



PA + Pm + 3H (p m + p m ) = . (2.8) 
An exchange of energy between the matter and the vacuum takes place. 

• Model III: G ^ and pa =const.: 

G( Pm + p A ) + G[p m + 3H( Pm + Pm )\ = . (2.9) 

Since G does not stay constant here, this equation implies a non-conservation of matter. It 
could be solved e.g. for G, if p m and p\ would be given by some non-conservation ansatz. 

• Model IV: G^O and p A / 0: 

There are many possibilities here. We consider the simplest one by assuming the standard 
local covariant conservation of matter-radiation, i.e Eq. (|2.7p . Eq. (j2.6|) leads to: 

( Pm + PA)G + Gp A = 0. (2.10) 

This situation is complementary to the previous one. Here the dynamical interplay is between 
G and p\, whereas p m is also time evolving, but decoupled from the feedback between G and 
PA- 

• Model V: Another possibility with G/0 and p\ ^ is the case that there is no matter in 
the universe: p m = 0. Then Eq. (j2.6[) implies G p\ =const. This does not exclude that both 
parameters can be time evolving while the product remains constant. This situation could 
only be of interest in the early universe, when matter still did not exist and only the vacuum 
energy was present. 

Only in the class of Models I and IV matter is covariantly self-conserved, i.e. matter evolves 
according to Eq. (|2.7p . In terms of the scale factor we find: 

p' m (a) + -(1 + cu m ) p m (a) = . (2.11) 

The prime indicates d/da. Its solution can be expressed as follows: 

p m (a) = p° m a-^ 1+ ^ = pi (1 + ,)3(i+^) . (2.12) 

We have expressed the result (|2.12p in terms of the scale factor a = a(t) and the cosmological 
redshift z = (1 — a) /a. 
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We shall focus on Models II, III and IV. Each of these models stands for a whole class of 
possible scenarios. One has to introduce more specifications before being able to perform concrete 
calculations. The variation of the "fundamental constants" (e.g. p\, G) could emerge as an 
effective description of some deeper dynamics associated to QFT in curved space-time, e.g. in 
quantum gravity or in string theory. This should provide definite time/redshift-evolution laws 
PA = Pa{z) , G = G(z). Examples will be discussed in the next sections. 

Other fundamental parameters could also be variable. The fine structure constant might change 
in time/redshift - see e.g. \26\ 127] . However positive evidences [28] are questioned [29]. The 
possibility that the fundamental QCD scale parameter Aqcd of the strong interactions could also 
be time-evolving (hence redshift dependent) is of special interest (see sections 0] and [5] for details). 
This could lead to the non-conservation of matter in the universe. In this paper we discuss the 
possibility that this non-conservation of matter might be related to the cosmological matter non- 
conservation. This would lead to a departure from the standard cosmological scenario. 



3 Running vacuum energy and the coupling of gravity 

The running of the vacuum energy and/or the gravitational coupling is expected in QFT in curved 
space-time [30^ [3"T] , see also [20] and references therein. Running couplings in flat QFT provide 
a useful theoretical tool to investigate theories as QED or QCD. Here the corresponding gauge 
coupling constants run with the typical energy of the process. 

In the universe we expect that the running of p\ and G is associated with the typical energy of 
the classical gravitational external field linked to the FLRW metric. Here the Hubble rate H will 
set the scale, since it is related to the non-trivial structure of the FLRW background. The universe 
in an accelerated expansion {H ^ 0, H ^ 0) is a space-time with dynamical intrinsic curvature: 

R = -s(°L + &\ = -\2H 2 -QH. (3.1) 
\ a a z ) 

In the effective action of QFT in curved space-time [32] p\ and G should be effective couplings 
depending on a mass scale p. This scale parameterizes the various quantum effects from the matter 
fields. In some cases the vacuum energy and the gravitational coupling can be represented as a 
power series of p. The rates of change are given by: 

A 2k p 2k = A + A 2 p 2 + A±p 4 + ..., (3.2) 

fc=0,l,2,... 

= Yl B 2k p 2k = B + B 2 p 2 + B±p A + .... (3.3) 

fc=0,l,2,... 

Such a "running" of p\ and G with p reflects the dependence of the leading quantum effects on 
a cosmological quantity £ associated with p, hence p\ = pa(0 and G = G(£). In cosmology we 
expect that the physical scale £ could be the Hubble rate H(t), or the scale factor a(t) [30], which in 
most of the cosmological past also maps out the evolution of the energy densities with H. We will 
concentrate here on the setting p = H, which naturally points to the non-trivial curvature of the 
background - Eq (I3.ip - and also to the typical energy of the FLRW "gravitons" attached to the 



i = £ 



d 

dlnp 2 \G(p) 
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quantum matter loops contributing to the running of pa and G~ l in a semi-classical description 
of gravity. The coefficients A 2 k , B 2 ^ receive contributions from boson and fermion matter fields 
of different masses Mj. The series (|3.2p becomes an expansion in powers of the small quantities 
H/Mi (see Eq. f|3.4[) below). Only even powers of H are involved, due to the general covariance 
of the effective action |3D1I5I]R These expansions converge very fast for p = H, since H/Mi 1 
for any ordinary particle mass. No other iJ 2n -terms beyond H 2 (not even H 4 ) can contribute 
significantly on the r.h.s. of equation (|3.2p at any stage of the cosmological history below the 
GUT scale M x < M P . We find: 



dp A (p) 1 



dlnfi 2 (47r) 



n 2 p 2 + 0(p A ). (3.4) 



We have omitted the Aq term - it would be of order Mf. This would produce a too fast running 
of p\. This can also be derived from the fact that all known particles satisfy p < Mi (for p = H). 
None of them is an active degree of freedom for the running of p\, and only the sub leading terms 
are available. Approximately we obtain a simple expression: 

p A (H) =n + n 2 H 2 . (3.5) 

In view of the boundary condition pa(Hq) = p A it is convenient to rewrite the coefficients of (|3.5p : 

n o = P°A-^M 2 H 2 , n 2 = ^-M 2 . (3.6) 

We have defined the dimensionless parameter 

v = Ly Ci ^. (3.7) 

6vr £-< l M 2 V ; 

i=f,b r 

The sum runs over fermions (/) and bosons (b) contributing to the loop. The parameter v provides 
the main coefficient of the one-loop /3-function for the running of the vacuum energy. The generic 
expression (j3.7[) adopts a concrete form with coefficients Cj, depending on the effective action of 
the underlying QFT (see e.g. [30]) . The parameter v can have any sign a = ±, depending on 
whether bosons or fermions dominate. 

It is convenient to write (|3.7p as follows: 

a M 2 . . 

v = T . 3.8 

6vr M 2 K ' 

Here M 2 = \ X^i=/ b c i^i \ ls an effective mass squared representing all the particles contributing 
to the running after counting their multiplicities. For M = Mp we have \v\ = 0(1O -2 ). In general 
we expect that the set of Mj includes masses of some GUT theory with a mass scale Mx ~ 10 16 
GeV (M ~ M x < M P ). A natural estimate is in the range v = 10 -6 - 1CT 3 [30] , 

If we would instead take the string scale as the characteristic GUT scale [341 135] . then M/Mp ~ 
10 -2 , and \u\ could move to the upper range 10~ 3 . For v = we have n 2 = in Eq. (|3. 5|) . In 



3 In practice, if one tries to fit the data with a time dependent CC term which is linear in the expansion rate, i.e. 
of the form A oc H , the results deviate significantly from the standard ACDM predictions [33]. 
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this case the vacuum energy remains strictly constant at all times: p A = p\ , and we recover the 
standard situacion of the ACDM model. For non- vanishing v the evolution law (|3. 51 leads to: 



p h{ H)=pl + ^-M 2 {H 2 -H 2 ). (3.9) 

The expansions (|3.2p - (|3,3p are correlated by the Bianchi identity (|2.6|) . If p = p(t) is a well defined 
invertible function, dp/dt 7^ - as it is in the case with p = H(t) - we must have 



dG dp a 

~T [Pm + PA) + G — h G 

dp dp 



dpm. 3 da 

-J h - [Pm + Pm) -7- 

dp a dp 



0. (3.10) 



This expression shows, that the dynamical dependence of pa and G may not be in the cosmic 
time t (as in many phenomenological models in the literature |24|). but in p. There is a possible 
connection of the evolution of p\ with the quantum effects of QFT in a curved background, i.e. 
with the running pa(p) in an expanding universe [20J. Since the quantum effects on G and p\ must 
satisfy the above differential constraint, they must be correlated. If we assume that p\ evolves as 
indicated in (j3.9|) . the corresponding running of G must fulfill ()3.10p . But this is still not enough 
to determine G = G(H) explicitly, since it depends on whether matter is conserved or not. Then 
one has to have a specific ansatz for the matter non-conservation equation. 

We consider two possibilities. We assume that matter is conserved, as in Model IV of the 
previous section. The term in brackets on Eq. (|3.10p vanishes - see ([2. lip . Using Friedmann's 
equation f|2.3j) and (|3.9p . we are left with: 

3 H 2 dG „ 3z/ 9 , 
87rG dH An 

G(H) = G ° - - „, . (3.12) 

V 1 1 + v\r(H 2 /H 2 ) V ; 

Here we have defined Go = 1/Mp, the current value of G, i.e. Go = G(H ). From (|3.12p we find: 

d 1 - r2 



After integration we obtain: 



i*B>a-" M *- (313) 

Thus (13.121) is the solution of (13. 3p . when we take only the leading term in the expansion, which does 
not depend on p = H. This is consistent, since 1/G is a large quantity and must be dominated by 
this term. The quantity p\, which in contrast is a much smaller quantity, can not be dominated by 
Aq ~ Mf, but rather by the next-to- leading term, which is proportional to H 2 . The leading term 
in each case dominates the corresponding running equation. Higher order corrections (involving 
more powers of H) are possible, but they are negligible in view of the current value of H. 

We mention another simple case, where matter is not conserved. We write dG/dp = G'(a) dp/ da, 
dp\/dp = p' A (a) dp/da, and dp m /dp = p' m (a) dp/da. Assuming that p = p(a) is a well defined 
invertible function (which is indeed the case, when p = H) we have dp/da 7^ 0. If G is constant, 
Eq. (|3.10p simplifies again: 

3 

p' A (a) + p' m (a) + -(1 + u m ) p m (a) = . (3.14) 
a 
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H 2 (a) 



(3.17) 



This result is consistent with (|2.8p . The running of the vacuum energy is due to the non- 
conservation of matter. The solution is well-known (see [20] and references therein). The cor- 
responding matter non-conservation law is: 

p m {a) = pl ia -^ 1+ ^ l -»K (3.15) 

The associated running of the vacuum energy density as a function of the scale factor is given by: 

p A (a) = pi + [ a -3(i+"™)(i-*>) _ jj . ( 3-16 ) 

The equations (|3.15p and (|3.16j) do satisfy (|3.14j) and the boundary conditions p m (a = 1) = p Q m 
and p\ (a = 1) = p A are fulfilled for the present universe. The running vacuum law (|3.16p is a 
consequence of the original equation (]3.9p . The consistency of these two formulae implies that the 
invertible function p = p{a) (i.e. H = H(aj) is given by: 

8itG 
3(1-!/) 

Here p® is the present value of the critical density. For v = we obtain a dilution law for the 
matter density of the form (I2.12p . i.e. ~ a -3 for non-relativistic and ~ a -4 for relativistic matter; 
also a constant p\ = p A , and the canonical form for the Hubble expansion rate H = H(a). 

These models are compatible with the observational data, both on the Hubble expansion (e.g. 
from SNIa+BAO) as well as on structure formation (power spectrum, growth factor and CMB) 
for values of the relevant parameter (|3.7p up to \u\ ~ 10~ 3 - see [331 ES] for details. We shall come 
back to this cosmological input in sect.O 

We note the coincidence of this order of magnitude estimate for v with its theoretical ex- 
pectations for being a f3- function coefficient of p\. The generalizations of these running vacuum 
models are possible at a similar level of phenomenological success, see e.g. [37]. As also shown 
in this reference, alternative dynamical models (such as the so-called entropic-force cosmologies) 
are not successful, although they have many elements in common. There exist time evolving vac- 
uum models, which can help to cure the old cosmological constant problem and the coincidence 
problem 0SI . 

Thus there exists an interesting class of cosmological models with time evolving vacuum energy 
which are phenomenologically acceptable, but not every phenomenological model can be success- 
fully tested (in this respect we have also mentioned the unsuccessful cosmologies with vacuum 
energy linear in H - see [33] and references therein) . 



4 Time evolving masses in the Standard Model of particle physics 

In this section we discuss experiments on the time variation of the fundamental constants of Nature. 
We suggest that they could be related to matter non-conservation. 

4.1 The Oklo phenomenon 

There are experiments which suggest that the fine structure constant a em has not remained con- 
stant throughout the cosmic evolution. There are many independent observations suggesting this 
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possibility [2HJ [27J . We also mention the "Oklo phenomenon" [T21 HH IS] • It is related to the 
natural fission reactor (the Oklo uranium mine) in Gabon (West Africa), first discovered in 1972 
by the French Commissariat a l'Energie Atomique. This natural reactor operated nearly 2 billion 
years ago for a period of some two hundred thousand years at a power of ~ 100 Kw. The data 
correspond to a process that occurred at the redshift z — 0.16 (for the typical values h ~ 0.70, 
£l° M ~ 0.27, ~ 0.73 of the cosmological parameters). This is the redshift at which we may be 
sensitive to variations of the fundamental constants. 

The fraction of 235 [7 in the Oklo site has decreased since then from 3.68% to 0.72%. This 
depletion with respect to the current standard value is a proof of the past existence of a spontaneous 
chain reaction. Water from river Oklo provided the moderator for the neutrons. One of the nuclear 
fission products is the Samarium's isotope 149 Smg2 which upon neutron capture becomes the excited 
isotope of the same element 150 Sni62: 

149 Sm 62 +n^ 150 Sm 62 + 7 . (4.1.1) 

The sustained fission chain at the Oklo mine leads to the process (|4.1.ip . The relatively light 
isotope 149 Sni62 is not a fission product of the 235 ?7, so the reaction (|4.1.ip took place in the 
natural ores of Oklo. It was observed that the ratio of isotopes 149 Sm62/ 147 Sni62 in samples of 
Samarium in these ores is 0.02, while in normal Samarium is 0.9. The depletion shows that the 
reaction (|4.1,ip took place for a long time in the Oklo reactor. 

The cross section of the neutron capture ()4.1.1|) depends on the energy of a resonance at 
E r = 97.3 meV and is well described by the Breit-Wigner formula: 

gnh 2 r w r 7 
a{E) -2^E(E-E r y + ryA- (4L2) 

Here g = 9/16 is a spin-dependent statistical factor, T is the total width, i.e. the sum of the 
neutron partial width (r n = 0.533 meV) and of the radiative partial width (r 7 = 60.5 meV). In 
order to estimate the cross-section in a more realistic way, one has to thermal average the above 
Breit-Wigner formula, using the geophysical conditions at the Oklo site. From here one can infer 
the uncertainty in the resonance energy, 5E r , which is set equal to _E° kl ° — E®, where E° kl ° is the 
value of the resonance during the Oklo phenomenon and E® is the possibly different value taken 
today. From the mass formula of heavy nuclei the change in resonance energy is related to a cm 
through the Coulomb energy contribution: 

SE r = -1.1 ^MeV. (4.1.3) 

From the estimates on 5E r (ranging from a dozen meV to a hundred MeV [141 115]) one infers from 
(|4.1.3|) a tight bound on the time variation of the fine structure constant of order d em /a cm ~ 
10~ yr" 1 . This is comparable to the best bounds from atomic clocks |26[ 127] . But the debate 
continues on the reliability of the data obtained in the Oklo mine. Even if the corresponding 
bound, obtained on the time variation of the electromagnetic coupling, is eventually validated, the 
Oklo phenomenon cannot easily provide information on the time variation of the strength of the 
nuclear interaction, since it is sensitive only to dimensionless ratios of nuclear quantities. It cannot 
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be used to extract a possible variation of the QCD scale parameter Aqcd- This is essential to 
establish a link between the time variation of fundamental nuclear and particle physics constants 
with the the corresponding variation of the vacuum energy density in the the cosmic expansion. 

4.2 Time variation of the fundamental QCD constant: implications for the 
nucleon mass and the nuclear masses in the universe 

It has been argued that the fundamental QCD scale parameter Aqcd could vary much faster than 
a em \10\ [TTj 112] . This change would be related to a corresponding change of the nucleon mass. 
Within the context of QCD the nucleon mass and the other hadronic masses are determined by 
the value of the QCD scale parameter Aqcd- The leading contribution to the nucleon mass can 
be expressed as tun — cqcdAqcd, where cqcd is a non-perturbative coefficient. The masses of 
the light quarks m u , m d and m s also contribute to the the proton mass, although by less than 
10% only. There is also a small contribution from electromagnetism. Let us take for instance the 
proton mass m p ~ 938 MeV. It can be computed from the QCD scale parameter Aqcd, the quarks 
masses and the electromagnetic contribution: 

m p = cqcdAqcd + c u m u + c d m d + c s m s + c cm AQ C D 

= (860 + 21 + 19 + 36 + 2) MeV. (4.2.1) 

The QCD scale parameter is related to the strong coupling constant a s = gi?/(47r). To lowest 
(1-loop) order one finds: 

1 47T 

a s (ji R ) = v = v , (4.2.2) 

A) ^ (Aq CD //4) (11 -2^/3) In (^/A QCD ) 

where [i R is the renormalization point and /3q = — bo = —(33 — 2ny)/(127r) (rif being the number 
of quark flavors) is the lowest order coefficient of the /3-function. 

The QCD scale parameter Aqcd has been measured: Aqcd = 217 + 25 MeV. When we embed 
QCD in the FLRW expanding background, the value of Aqcd need not remain rigid anymore. 
The value of Aqcd could change with H, and this would mean a change in the cosmic time. If 
Aqcd = Aqcd(-£0 is a function of H, the coupling constant a s = a s (/j, R ;H) is also a function of 
H (apart from a function of fi R ). The relative cosmic variations of the two QCD quantities are 
related (at one-loop) by: 

1 da s {fx R \H) _ 1 
a s dH In (ji R / Aqcd) 

If the QCD coupling constant a s or the QCD scale parameter Aqcd undergo a small cosmological 
time shift, the nucleon mass and the masses of the atomic nuclei would also change in proportion 
to Aqcd- 

The cosmic dependence of the strong coupling a s ((i R ;H) can be generalized to the other 
couplings Qj = 0-i{fJL R \ H) [11] , In a grand unified theory these couplings converge at the unification 
point. Let da,i be the cosmic variation of a, with H. Each of the oti is a function of fi R , but the 



Aqcd 



dH 



(4.2.3) 
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expression a i 1 {doti/on) is independent of /iR. One can show that the running of a cm is related to 
the corresponding cosmic running of Aqcd as follows: 

1 dA QCD (HY 



a 



1 da cm (fi R ; H) _ 8 a em (fj, R ;H)/a s (fi R ;H) 
~ 3 



dH 



In (mr/Aq CD ) 



A 



QCD 



dH 



(4.2.4) 



At the renormalization point fj> R = Mz, where both a cm and a s are well-known, one finds: 



1 da cm (^ R ;H) 
a Pm dH 



0.03 



dA QCB (H) 



A 



QCD 



dH 



(4.2.5) 



Thus the electromagnetic fine structure constant runs more than 30 times slower with the 
cosmic expansion than Aqcd- Searching for a cosmic evolution of Aqcd is much easier than 
searching for the time variation of a em . 



4.3 Time evolution of the proton - electron mass ratio 

We consider the mass ratio: 



P'pe — 



(4.3.1) 



This ratio is known with high accuracy: /i pe = 1836.15267247(80) |39j . Since a change of Aqcd 
would not affect the electron mass, the mass ratio (|4.3.ip would change during the cosmological 
evolution. 

First we consider astrophysical tests. The spectrum of H2 provides a direct operational handle 
to test possible variations of (|4.3.ip . Particularly significant is the study of Ref.[19], based on 
comparing the H2 spectral Lyman and Werner lines, observed in the Q 0347-383 and Q 0405-443 
quasar absorption systems, with the laboratory measurements. 

The result indicates, that n pe could have decreased in the past 12 Gyr, corresponding to a 
relative time variation of 

^ = (-2.16 ± 0.52) x 10" 15 yr- 1 . (4.3.2) 

It has been pointed more recently by other authors [30] that this measurement may suffer from 
spectral wavelength calibration uncertainties, and the reanalysis of the time variation would show 
a significance at the 1 a level only. 

Now we consider laboratory tests, using atomic clocks. According to our estimate (|4.2.5p . the 
largest effect is expected to be a cosmological redshift (hence time variation) of the nucleon mass, 
which can be observed by monitoring molecular frequencies. These are precise experiments in 
quantum optics, e.g. obtained by comparing a cesium clock with 1S-2S hydrogen transitions. In 
a cesium clock the time is measured by using a hyperfine transition^. Since the frequency of the 
clock depends on the magnetic moment of the cesium nucleus, a possible variation of the latter is 
proportional to a possible variation of Aqcd- A hyperfine splitting is a function of Z a em (Z being 
the atomic number) and is proportional to Z a^d^N / /j,B)(m e /m p ) Roo, where Roo is the Rydberg 

4 Recall that the cesium hyperfine clock provides the modern definition of time. In SI units, the second is defined 
to be the duration of 9.192631770 x 10 9 periods of the transition between the two hyperfine levels of the ground 
state of the 133 Cs atom 
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constant, fj,pf is the nuclear magnetic moment and ps = eh/2m p c is the nuclear magneton. We 
have ftN/fJ'N tx — Aqcd /Aqcd- The hydrogen transitions are only dependent on the electron 
mass, which we assume to be constant. The comparison over a period of time between the cesium 
clock with hydrogen transitions provides an atomic laboratory measurement of the ratio (|4.3.ip . 
The most recent atomic clock experiment at the MPQ (Max-Planck-Institut fur Quantenoptik) at 
Garching near Munich gives a limit |16| : 

<10 _14 yr _1 . (4.3.3) 

Since the proton mass is given essentially by Aqcd , as indicated by Eq. (|4.2.ip , we have m p ~ 
ca qcd Aqcd- The corresponding time variation of the ratio (I4.3.ip would be: 

Ppe 

Thus the atomic clock result (|4.3.3p would indicate a time variation of the ratio p pc , which is 
consistent (in absolute value) with the astrophysical measurement (|4.3.2p . The result above implies 
also a bound for a possible time variation of the light quark masses: 

<10- 14 yr _1 - (4-3.5) 

5 Dynamical dark energy and a cosmic link with nuclear and 
particle physics 

The time evolution of the fundamental "constants" p\ and G of gravity could be related to the 
time variation of the fundamental "constants" in nuclear and particle physics. In some models one 
can have matter conservation even though p\ is running, but at the expense of having a running 
G as well - confer Model IV of sect. [2] and Eq. ()2.10p . In an alternative class of models, G runs 
thanks to the non-conservation of matter, as in Model III of sect. but then p\ stays fixed. If G 
stays fixed and p\ is evolving, there is a transfer of energy from matter into the vacuum, or vice 
versa - cf. the Model II class of sect. [2] and Eq. (|2.8p . The various classes of cosmological scenarios 
are interesting, but the last two could help us to understand the potential cosmic time variation of 
the fundamental "constants" of nuclear and particle physics, such as the QCD scale, the nucleon 
mass and the masses of nuclei. 

5.1 Non-conservation of matter at fixed G 

First we consider the class of scenarios denoted as Model II. Let p° M be the total matter density of 
the present universe, which is essentially non-relativistic (oj m ~ 0). The corresponding normalized 
density is = p® M /p° c — 0.27, where p® is the current critical density. Similarly, = ~ 
0.73 is the current normalized vacuum energy density, for flat space. If p\ evolves with the Hubble 
rate in the form indicated in Eq. (|3.9p . the non-relativistic matter density and vacuum energy 



Aqcd 
Aqcd 



m p 




Aqcd 






Aqcd 



< 10- 14 yr _1 



(4.3.4) 



m q 
m a 
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density evolve with the scale factor, given in (|3.15p and (|3.16p . Expressing the result in terms of 
the cosmological redshift z = (1 — a)/a, we find: 

p M (z;u)=p° M (l + z) 3 ^, (5.1.1) 

and 

PA(z)=pi + ^ [(1 + ~ l] ■ (5.1.2) 

The crucial parameter is v, which we have introduced in sect.[3l It is responsible for the time 
evolution of the vacuum energy. From Eq. (|5.1.1[) we confirm, that it accounts also for the non- 
conservation of matter, since it leads to the exact local covariant conservation law (|2.11|) . For 
non-relativistic matter we find: 

Pm{z) = p° M (1 + zf . (5.1.3) 

5pM = Pm{z\v) — Pm(z) is the net amount of non-conservation of matter per unit volume at a 
given redshift. This expression must be proportional to u, since we subtract the conserved part. 
At this order we have bpu = — 3v p° M {l + zf \n(l + z). We differentiate it with respect to time and 
expand in u, and finally divide the final result by pm- This provides the relative time variation: 

^ = 3v (l + 31n(l + z)) H + 0{u 2 ). (5.1.4) 
Pm 

Here we have used z = (dz/da)a = (dz/da)aH = —(1 + z)H. Assuming relatively small values of 
the redshift, we may neglect the log term and are left with: 

Spm 

PM 

From (|5.1.2p we find: 



3vH. (5.1.5) 



— ~ —2>u (1 + zf H + 0(u 2 ) • (5-1.6) 
PA ft A 

It is of the same order of magnitude as (|5.1.5|) and has the opposite sign. Let us compare the 

theoretical expression (|5.1.5p with the experimental results f|4.3.2[) and (|4.3.4p . described in the 

previous section. Taking the current value of the Hubble parameter as a reference, Hq = 1.0227 hx 
10 -io yr -i f where h ^ q 7Q ) 

we obtain \v\ < O(10 4 ) for the most conservative case. It is a rather 
tight bound, in accordance with the QFT expectations in sect.[3J 

What is the role played by the running vacuum energy (|5.1.2p ? Its evolution in combination 
with the non-conservation of matter affects many relevant cosmological observables, which have 
currently been measured with high precision. From a detailed analysis of the combined data on 
type la supernovae, the Cosmic Microwave Background (CMB), the Baryonic Acoustic Oscillations 
(BAO) and the structure formation data a direct cosmological bound on v has been obtained in 
the literature [331 ES] : 

l^cosm. < O ( 10 -3) j (Model II sect. [2} . (5.1.7) 

It is consistent with the theoretical expectations. In the next section we analyze another model 
which can also accommodate matter non-conservation in the form (|5.1.1|) . but at the expense of a 
time varying G. We compare it with a similar model, where matter is conserved. 
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5.2 Non-conservation of matter at fixed p A 



Within the class of scenarios indicated as Model III of sect. [3] the parameter p\ remains constant 
(PA = Pa) an< ^ G * s var i a ble. This is possible due to the presence of the non self-conserved matter 
density (|5.1.ip . Trading the time variable by the scale factor, we can rewrite Eq. (|2.9p as follows: 



G'{a) [p M (a) + p° A ] +G(a) 



Pm(«) + -Pm{o] 



0. 



(5.2.1) 



The primes indicate differentiation with respect to the scale factor. We insert equation (]5.1.ip in 
(|5.2.ip . integrate the resulting differential equation for G(a) and express the final result in terms 
of the redshift: 



GO) = Go 



u/(l-v) 



(5.2.2) 



Here Go = 1/Mj, is the current value of the gravitational coupling. The previous equation is 
correctly normalized: G(z = 0) = Go, due to the cosmic sum rule in flat space: + = 1. For 
v = the gravitational coupling G remains constant: G = Gq. Since p\ is constant in the current 
scenario, the small variation of G is entirely due to the non-vanishing value of the z^-parameter in 
the matter non-conservation law (|5.1.ip . This leads to the dynamical feedback of G with matter[§. 
For the present model Friedmann's equation (|2.3|) becomes: 



H\z) 



8ttG(z) 



P°m (1 + 



Hi 



G(z) 

Gn 



n° M (i + zf^ + 



Combining (|5.2.2p and (|5.2.3p . we find the Hubble function of this model in terms of z: 

1/(1-") 



H\z) 



Hi 



O 



and we obtain: 



G(z) 
Go 



fl +2 )3(l-") +n 



H\z) 
Hi 



(5.2.3) 



(5.2.4) 



(5.2.5) 



Since v is presumably small in absolute value (as in the previous section), we can expand (|5.2.5p 
in this parameter: 

G{H) ~ G ^1 + u In ^ + 0{v 2 )^j . (5.2.6) 

At leading order in v this expression for the variation of G is identical to the one found for Model 
IV of sect. [21 see Eq. (|3.12p . except for the sign of v. The equation (|5.2.6p allows us to estimate 
the value of the parameter v by confronting the model with the experimental data on the time 
variation of G. Differentiating (15.2.6P with respect to the cosmic time, we find in leading order in 
v. 



G 
G 



2u' 



H 



-2{l + q)vH, 



(5.2.7) 



H 

where we have used the relation H = — (1 + q)H 2 , in which q = —d/(aH 2 ) is the deceleration 
parameter. From the known data on the relative time variation of G the bounds indicate that 



5 This feedback can also be conceived in the context of gravitation holography [41] if one also takes as a starting 
point the matter non-conservation law (J57T7TJ) . This law was first suggested and analyzed in [42] and later on in [43] . 
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\G/G\ < 10 12 yr 1 [261 127). If we take the present value of the deceleration parameter, we have 
q = 3tt° M /2 - 1 = -0.595 ~ -0.6 for a flat universe with Q° M = 0.27. It follows: 



G 



G 



< 0.8\v\H. (5.2. 



Taking the current value of the Hubble parameter: Hq ~ 7 x 10~ n yr -1 (for h ~ 0.70), we obtain 
\v\ < 10 -2 . The real value of |^| can be smaller, but to compare the upper bound that we have 
obtained with observations makes sense in view of the usual interpretation of v in sect. [3] and the 
theoretical estimates indicated there. The constraints from Big Bang nucleosynthesis (BBN) for 
the time variation of G are more stringent and lead to the improved bound: 

M BBN < 10" 3 , (Model III sect. ED . (5.2.9) 

This bound can be obtained by adapting the study of Ref. [44] , which was made for Model IV 
of sect. El Since Models III and IV share a similar kind of running law for the gravitational 
coupling (except for the sign of v) - confer equations ()3. 12[) and ()5.2.6|) — , we can extract the 
same bound for \v\ in the two models following the method of sect. 5.2 of Ref. [33] and references 
therein, particularly [45] . The final result is Eq. (|5.2.9|) . The cosmological data from different 
sources furnish about the same upper bound on \v\ for the two running models where matter is 
non-conserved, i.e. Models II and III of sect. El In both cases the upper bound on \v\ is ~ 10~ 3 , 
as shown by equations (|5.1.7p and (|5.2.9p . 

The previous bounds on \v\ for Models II and III are completely general (meaning that they 
apply to all forms of matter), since they are obtained from cosmological data tracing the possible 
evolution of p\ and G, respectively. But these cosmological bounds are weaker than those that 
follow, if we interpret v as a matter non-conservation parameter. Since matter is indeed non- 
conserved in both of these models, Eq. (|5.1.5|) and the lab bound (|4.3.3|) do apply in the present 
case, but only if the non-conserved matter is of nuclear nature. In this case we obtain the stronger 
constraint 

l^iab. < O ( 10 -4) , (Models II and III sect. ED . (5.2.10) 

But if the non-conserved matter is dark matter, then only the weaker (purely cosmological) bound 
(|5.2.9p is valid (see the next section for a detailed discussion on the distinct contributions from 
nuclear matter and dark matter). 

Despite |G| varies with time in a comparable way in Models III and IV, the stronger bound 
(|5.2.10p does not apply for Model IV, since matter is conserved in it and hence Eq. ()5.1.ip does not 
hold for this model. Only the pure BBN cosmological bound (|5.2.9p is applicable in this case. This 
primordial nucleosynthesis bound on Model IV coincides with an independent bound obtained for 
this model from type la supernovae, the Cosmic Microwave Background, the Baryonic Acoustic 
Oscillations and the structure formation data (cf. [36] for details). For Model IV two independent 
cosmological bounds (BBN plus the current cosmological data) converge to the same result: 

j^BBN+cosm. < O ( 10 -3) j (Model IV sect. ED . (5.2.11) 
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Although the order of magnitude of the bounds on \ v\ are sometimes coincident for different models, 
they are different. For example, Model IV cannot - in contrast to Models II and III - be used to 
explain the possible time variation of the fundamental constants of the strong interactions and the 
particle masses. It can only be used to explain the time variation of the cosmological parameters 
PA and G in a way which is independent from the microphysical phenomena in particle physics 
and nuclear physics. 

Finally, we note that the above cosmic changes in the values of the proton to electron mass 
ratio and G or pa can be written in terms of dimensionless quantities (in natural units). For 
example, for Model II (where G is fixed and p\ is variable) we can define the dimensionless 
quantity A = A./m p = 8tt G p\/m p . Then, 

1 dX pa ih„ 

-— = — -2-2- oc uH, (5.2.12) 
A at pa m p 

because both terms on the r.h.s. are proportional to v (cf. sect. 15.1ft . Similarly, for Model III 
(where p\ is fixed and G is variable) we can construct the dimensionless quantity GmJ, Its 
relative variation is also proportional to u: 



1 d(Gm) G 



m 



+ 2^-otvH. (5.2.13) 



Grrip dt G m, 

5.3 Non-conservation of baryonic matter versus dark matter and the cosmic 
evolution of Aqcd 

Here we focus on the impact of the cosmological Models II and III of sect.[2]on the non-conservation 
of matter in the universe. In the previous section we have considered bounds on the "leakage 
parameter" v within the class of these models based on the non-conservation matter density law 
(]5.1.ip . We must be careful in interpreting such a non-conservation law. For example, if we take 
the baryonic density in the universe, which is essentially the mass density of protons, we can write 
Pm = n p m pi where n p is the number density of protons and m® = 938.272013(23) MeV is the 
current proton mass. If this mass density is non-conserved, either n p does not exactly follow the 
normal dilution law with the expansion, i.e. n p ~ a~ 3 = (1 + z) 3 , but the anomalous law: 

n p {z) = rip (1 + z) 3< - 1_l/ ^ (at fixed proton mass m p = m p ) , (5.3.1) 

and/or the proton mass m p does not stay constant with time and redshifts with the cosmic evolu- 
tion: 

m p (z) = (1 + z)~ 3u (with normal dilution n p (z) = n° (1 + zf) . (5.3.2) 

In all cases it is assumed that the vacuum absorbs the difference (i.e. pa = P\( z ) "runs with the 
expansion" ) . The first possibility implies that during the expansion a certain number of particles 
(protons in this case) are lost into the vacuum (if v < 0; or ejected from it, if v > 0), whereas 
in the second case the number of particles is strictly conserved. The number density follows the 
normal dilution law with the expansion, but the mass of each particle slightly changes (decreases 
for v < 0, or increases for v > 0) with the cosmic evolution. 
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Here we adopt the second point of view, i.e. Eq. (|5.3.2p . We can interpret the tight bounds from 
the laboratory and cosmological observations summarized in sect. 0] as direct bounds on the cosmic 
time evolution of Aqcd (hence on m p and on the nuclei in the universe) . Since the contribution of 
the quark masses m u ,md and m s to the proton mass is small - cf. Eq. (|4.2.ip - we can approximate 
the proton mass by m p ~ cqcd Aqcd • It will be sufficient to take into account the leading effects 
of the time variation of m p through the corresponding effects in Aqcd • 

Since the matter content of the universe is dominated by the dark matter (DM), we cannot 
exclude that it also varies with cosmic time. Let us denote the mass of the dominant DM particle 
mx, and let px and nx be its mass density and number density, respectively. The overall matter 
density of the universe can be written as follows: 

PM = PB + PL + PR + PX = (n p m p + n n m n ) +n e m e + PR + n x m x 

~ n p m p + n n m n + n x mx • (5.3.3) 

Here n p ,n n ,n e ,nx (m p ,m n ,m e ,mx) are the number densities (and masses) of protons, neutrons, 
electrons and DM particles. The baryonic and leptonic parts are ps = n p m p + n n m n and pi = 
n e m e respectively. The small ratio m e /m p ~ 5 x 10~ 4 implies that the leptonic contribution to 
the total mass density is negligible: pi <C pb- We have also neglected the relativistic component 
Pr (photons and neutrinos). 

If we assume that the mass change through the cosmic evolution is due to the time change 
of m p , m n and mx, we can compute the mass density anomaly per unit time, i.e. the deficit or 
surplus with respect to the conservation law, by differentiating (|5.3.3|) with respect to time and 
subtracting the ordinary (i.e. fixed mass) time dilution of the number densities. The result is: 

SpM = n p m p + n n rh n + nx mx ■ (5.3.4) 

The relative time variation of the mass density anomaly can be estimated as follows: 



5pM _ n p m p + n n m n + n x m x ^ n p m p + n n m n + n x m x ( _ n p m p + n n m Ti 
p M n p m p + n p m p + n x m x n x m x \ n x m x 



(5.3.5) 



The current normalized DM density ^£) M = Px/Pc — 0.23 is significantly larger than the cor- 
responding normalized baryon density Q B = Pb/Pc — 0.04. Therefore nxmx is larger than 
n p m p + n n m n by the same amount. If we assume m n = m p , we find approximately: 

5p M n p m p ( n n Q B \ m x ( ^b \ „ 

1 + — " o + — 1 - o • ( 5 - 3 - 6 ) 



PM n x m x V n P ^dm/ m x \ ^dm. 
In the approximation m n = m p we can rewrite the prefactor on the r.h.s of Eq. (|5.3.6j) as follows: 



&B rh p ( n n /n p \ ^ tt B m p 



1-^). (5.3.7) 



n x m x ^dm m p \ l + n n /n p J 12dm m p \ n 

The ratio n n /n p is of order 10% after the primordial nucleosynthesis. Since Qb/^dm is also of 
order 10%, we can neglect the product of this term with n n /n p . When we insert the previous 



18 



equation into (|5.3.6[) . the two n n /n p contributions cancel each other. The expression 1 — J^b/^dm 
factorizes in the two terms on the r.h.s of Eq. (|5.3.6j) . The final result is: 

A M B \ " 1 SpM _ Mb rh p | ra x _ ^B A QC p | m x (5 3 8) 

V ^dm) Pm ^dm ni p ra x ^dm Aq C d mx 

We have used ra p ~ cqcdAqcd, the latter being accurate up to 10% corrections at most - see 

(|4,2,ip . Equation (|5.3.8p should be a good approximation (at most 10% corrections). 

The expression Spm/pm in Eq. (|5.3.8|) must be the same as the one we have computed in (|5.1.4p . 

if we consider the models based on the generic matter non-conservation law ()5.1.1[) . Therefore the 

two expressions should be equal, and we obtain approximately: 

o tt Aqcd t ra x „ n , 

3u eB H = 1 , (5.3.9) 

S'DM A QCD rax 

where we have defined 

^eff= 1 U /Q ■ (5.3.10) 
1 — ssb/' 'DM 

We have v % g ~ 1.2 zy. The differential equation (|5.3.9|) describes approximately the connection 
between the matter non-conservation law (I5.1.ip . the evolution of the vacuum energy density p\ 
(and/or G) and the time variation of the nuclear and particle physics quantities. Even if the DM 
does not change with the cosmic expansion, it is necessary to include it as a part of the total 
energy density of the universe. 

We assume that the dark matter particles do not vary with time, i.e. rax = 0, and only the 
cosmic evolution of Aqcd accounts for the non-conservation of matter. Trading the cosmic time 
for the scale factor through Aqcd = (i^Aqcd/^) o,H and integrating the resulting equation, we 
can express the final result in terms of the redshift: 

Aqcd(z) = Aqcd (1 + z)' 3 ^^)^ . (5311) 

For the protons we obtain: 

ra p {z) = m° p (1 + z)- 3 ^)^ _ (5.3.12) 

Here Aq CD and m® are the QCD scale and proton mass at present (z = 0). f2 DM and Q° B are the 
current values of these cosmological parameters. 

The presence of the factor 0^ / ^ DM in the power law makes eq. (|5.3.12p more realistic than eq. 
(|5.3.2p . In the case v = the QCD scale and the proton mass would not vary with the expansion 
of the universe, but for non- vanishing v it describes the cosmic running of Aqcd = Aqcd(^) and 
ra p = ra p (z). For v > (v < 0) the QCD scale and proton mass decrease (increase) with the 
redshift. This is consistent, since for v > {y < 0) the vacuum energy density is increasing 
(decreasing) with the redshift - cf. Eq. (I5.1.2P -, and it is smaller (larger) now than in the past. 

We can write down the variation of the QCD scale in terms of the Hubble rate H. With the 
help of Eq. (|3.17p equation (|5.3.1ip can be turned into an expression for Aqcd given explicitly in 
terms of the primary cosmic variable H : 



Aqcd(#) = Aqcd 



o 



1 - v H 2 tt° A 



O H 2 O 



(5.3.13) 
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with = Q B + v and v e g are involved in ([5.3. 13 j) . since they come from different sources. 

This equation satisfies the normalization condition Aqcd(-£^o) = Aqcd due to the cosmic sum rule 
for flat space: £79, + fi^ = 1. Using the previous equations and Eq. (|4.2.2p . we can obtain the 
corresponding evolution of the strong coupling constant a s with the redshift and the Hubble rate, 
i.e. a s (n R ;z) and a s (fi R ;H): 

1 1 O 

— r = —, ~r + 6 b -2M UcS l n (1 + z) . (5 . 3 .14) 

a s (/j, R ;z) a s (/i R ;0) il B 

Here a s (n R ; 0) is the value of a s (/j, R ; z) today (z = 0). Since bo > (cf. sect. l4.2h . we observe that 
for > (v < 0) the strong interaction a s (/i R ; z) decreases (increases) with z, i.e. with the cosmic 
evolution. We also find|§: 



1 1 +2feo ^|M^ ln 



1 - v H 2 



0° H 2 O 



(5.3.15) 



a s (fi R ;H) a s (n R ;Ho) Q° B 1 

Here a s (fi R ; Ho) is the current value of a s (ii R ; H). 

Above we have determined the strong coupling as a function of two running scales: one is the 
ordinary QCD running scale [i R , the other is the cosmic scale defined by the Hubble rate H, which 
has dimension of energy in natural units. The second term on the r.h.s. depends on the product 
of the two /5-function coefficients, the one for the ordinary QCD running (bo) and the one for the 
cosmic running (y oc f e ff). 

We find: 

i) for v = there is no cosmic running of the strong interaction, 

ii) for v > the strong coupling a s (fi R ;H) is "doubly asymptotically free". It decreases for large 

ji R and also for large H, whereas for v < the cosmic evolution drives the running of a s 
opposite to the normal QCD running, 

iii) the velocity of the two runnings is very different, because H is slowly varying with time and 

\u\ <S 1 and \u\ <S bo < 1. The cosmic running only operates in the cosmic history and is 
weighed with a very small /3-function. But it may soon be measured in the experiments with 
atomic clocks and through astrophysical observations. 

The previous equations describe not only the leading cosmic evolution of the QCD scale and 
the proton mass with the redshift and the expansion rate H of the universe, but they can account 
for the redshift evolution of the nuclear masses. For the neutron we can write approximately: 
m n ~ cqcd Aqcd • For an atomic nucleus of current mass Ma and atomic number A we have 
Ma = Z m p + (A — Z) m n — Ba, where Z is the number of protons and A — Z the number of 
neutrons, and Ba is the binding energy. Although Ba may also change with the cosmic evolution, 
the shift should be less significant, since at leading order the binding energy relies on pion exchange 
among the nucleons. The pion mass has a softer dependence on Aqcd: Tn n ~ \fm q Aqcd> due to 
the chiral symmetry. 



6 It is interesting to note that a similar running of a s with the cosmic expansion was pointed out in a different 
context by J.D. Bjorken in [46] . 
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In the previous approximations we have neglected the light quark masses m q . We can assume 
that the binding energy has a negligible cosmic shift as compared to the masses of the nucleons. 
In the limit where we neglect the proton-neutron mass difference and assume a common nucleon 
mass m® N at present, the corresponding mass of the atomic nucleus at redshift z is given at leading 
order by: 

M A {z) ~ Am% (1 + z)" 3 ^/^)^ _ Ba ■ (5.3.16) 

Although the chemical elements redshift their masses, a disappearance or overproduction of nuclear 
mass (depending on the sign of u) is compensated by a running of the vacumm energy p\, which 
is of opposite in sign, see (|5.1.6|) . 

Above we have described a simplified case, in which the nuclear matter evolves with the cosmic 
evolution as a result of the evolution of the fundamental QCD scale. In this scenario the light 
quark masses are neglected, and the DM does not participate in the cosmic time evolution. 

Alternatively we can assume that the nuclear matter does not vary with time, i.e. Aqcd = 0, 
and only the DM particles account for the non-conservation of matter. In general we expect a 
mixed situation, in which the temporal rates of change for nuclear matter and for DM particles 
are different: 

— = 3 z/qcd =3v x H. (5.3.17) 

Aqcd rn x 

We have defined the QCD time variation index, which is characteristic of the redshift rate of the 

QCD scale, while v x is the corresponding one for the DM. In this more general case we find: 

Aqcd W = A° CD (1 + z y****> mx{z) = m° x (1 + z)~ 3 vx . (5.3.18) 

We introduce the effective baryonic redshift index z/73: 

vb = ^qcd • (5.3.19) 

"DM 

The equations (|5.3. 181) satisfy the relation (|5.3.9I) . provided the coefficients ub and v x are related 
by 

v c s = v B + v x ■ (5.3.20) 

z/qcd is the intrisic cosmic rate of variation of the strongly interacting particles. The effective 
index vb weighs the redshift rate of these particles taking into account their relative abundance 
with respect to the DM particles. Even if the intrinsic cosmic rate of variation of Aqcd would be 
similar to the DM index (i.e. if z/qcd ^ v x)i the baryonic index (|5,3.19p would still be suppressed 
with respect to u x , because the total amount of baryon matter in the universe is much smaller 
than the total amount of DM. 

In this mixed scenario the mass redshift of the dark matter particles follows a similar law as 
in the case of protons (|5.3.12p . except now we have v c & — > ub- The proton would have the index 
z/qcd characteristic of the free (and bound) stable strongly interacting matter: 

m p (z) = m° p (1 + s)- 3 ("dm/Ob) *b = m (1 + z y^Qou _ (5.3.21) 

The DM particles have another independent index v x . The sum (|5.3.2U|) must reproduce the 
original index z/ c g oc u, which we associated with the non-conservation of matter. 
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Finally we consider the possible quantitative contribution to the matter density anomaly from 
the dark matter. The global mass defect (or surplus) is regulated by the value of the v parameter, 
but the contribution of each part (baryonic matter and DM) depends on the values of the individual 
components z/g and vx- We can obtain a numerical estimate of these parameters by setting the 
expression (|5,3.8|) equal to (|5.1.5|) . The latter refers to the time variation of the matter density 
PM without tracking the particular way in which the cosmic evolution can generate an anomaly 
in the matter conservation. The former does assume that this anomaly is entirely due to a cosmic 
shift in the masses of the stable particles. Taking the absolute values, we obtain: 

3\u cS \ H ~ 

Here we have used the experimental bound (|4.3.3p on the time variation of Aqcd- 

Several cases can be considered, depending on the relation between the intrinsic cosmic rates of 
variation of the strongly interacting particles and DM particles, ^qcd and vx- Since these indices 
can have either sign, we shall compare their absolute values: 

• 1) \v x \ < \vb\- 

This condition implies \v x \ <C I^qcdI- By demanding the stronger condition \vx\ <C \ v b\, 
we insure that the intrinsic QCD cosmic rate |^qcd| is much larger than the corresponding 
DM rate \vx\- We can neglect the rhx/mx term on the r.h.s. of ()5.3.22p . and we recover 
the equations (I5.3.11I) - (I5.3.15I) with u c g ~ vb- Using Hq ~ 7 x 10" 11 yr" 1 , we find: 

kxl^O, li/effl ~ \u B \ < 10" 5 , |i/qcd| < 5 x 10" 5 . (5.3.23) 

The bound on vb — v c s that we have obtained above can be compared with (|5.2.10p . The 
former (which is more stringent) is more realistic than the latter because here we have taken 
into account explicitly the suppression factor Q. B / ^dm of baryonic matter versus dark matter 
- and also the (small) difference between v and f e g. 

• 2) \vx\ ^ \vb\- 

Here we still have \vx\ smaller than |^qcd|> but the requirement is weaker. It follows: 
jfcff | - 2|i/x| - 2\v B \ = 2(O s /fi DM ) |^qcd|, and we find 

|z/ eff | < 2 x 10~ 5 , \v x \ ^ Wb\ < 10" 5 , Iz^qcdI < 5 x 10" 5 . (5.3.24) 

• 3) \u x \ ^ |^qcd|: 

The two intrinsic cosmic rates for strongly interacting and DM particles are similar, i.e. 
Aqcd / Aqcd and rh x / mx do not differ significantly. In this case Eq. (|5.3.22p leads to 

3\u eS \ H < f A + lj x 1(T 14 yr" 1 . (5.3.25) 

There are two sign possibilities (^qcd — ^ v x)i and we take the absolute value: 

l^cffl < ( 7p^- + 1 J |^QCd| ^ I^qcdI • (5.3.26) 

V'DM / 



4 A 



23 A 



qcd rn x 



QCD m x 



<23 Xl ° 



■ 14 yr-' + 



m x 
m x 



(5.3.22) 
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1 y 1 cosm 


M lab = Wb\ 


1 , , 1 cosm 

\ v x\ 


Model II 


1(T 3 (SNIa+BAO+CMB) 


1CP 5 (Atomic clocks+Astrophys.) 


10~ 3 


Model III 


1(T 3 (BBN) 


1CP 5 (Atomic clocks+Astrophys.) 


10 -3 


Model IV 


10 -3 (SNIa+BAO+CMB)+BBN 









Table 1: Upper bounds on the running index \v\ for the various models defined in sect.[2j Only 
for Models II and III a non- vanishing value of \u\ is related to non-conservation of matter and a 
corresponding time evolution of p\ and G, respectively. For these models, a part of v (viz. vb) is 
accessible to lab experiments, whereas the DM contribution (yx) can only be bound indirectly from 
cosmological observations (same cosmological bound as for the overall v). For Model IV matter is 
conserved, and a non- vanishing value of \v\ (only accessible from pure cosmological observations) 
is associated to a simultaneous time evolution of p\ and G - with no microphysical implications. 



We find: 

M < |"QCd| ^ Wx\ < 5 x 10" 5 . (5.3.27) 

4) |^qcd| < Wx\- 

Here the nuclear part is frozen. The non-conservation of matter is entirely due to the time 
variation of the DM particles. Eq. (|5.3.22|) gives: 

3uH^ (l-^-) • (5.3.28) 
m x V "DM / 

We have written this expression directly in terms of the original v parameter. In this case 
we cannot get information from any laboratory experiment on rhx/mx, but we do have 
independent experimental information on the original v value (irrespective of the particular 
contributions form the nuclear and DM components). It comes from the cosmological data 
on type la supernovae, BAO, CMB and structure formation. The analysis of this data |33|, I36| 
leads to the bound (|5.1.7p . which applies to all models, in which matter follows the generic 
non-conservation law (|5.1.ip and the running vacuum law (|3,9p — or the same matter non- 
conservation law and the running gravitational coupling law ()5.2.6p . as shown in Eq. (|5.2.9p . 
Since it depends on the cosmological effects from all forms of matter, it applies to the DM 
particles in particular. We find: 

|„ x [«»n<10- 3 . (5.3.29) 

This bound is significantly weaker than any of the bounds found for the previous scenarios 
in which the nuclear matter participated of the cosmic time variation. It cannot be excluded 
that the matter non-conservation and corresponding running of the vacuum energy in the 
universe is mainly caused by the general redshift of the DM particles. In this case only 
cosmological experiments could be used to check this possibility. If the nuclear matter also 
participates in a significant way, it could be analyzed with the help of experiments in the 
laboratory. For a summary of the bounds, see Table 1. 
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If in the future we could obtain a tight cosmological bound on the effective f e ff-parameter 
(|5,3,20p . using the astrophysical data, and an accurate laboratory (and/or astrophysical) bound on 
the baryonic matter part vb, we could compare them and derive the value of the DM component 
vx- If ^cfr an d vb would be about equal, we should conclude that the DM particles do not 
appreciably shift their masses with the cosmic evolution, or that they do not exist. If, in contrast, 
the fractional difference | (f e ff — ^sV^eS I would be significant, the DM particles should exist to 
compensate for it. 

6 Conclusions 

In this paper we have described theoretical models based on the assumption that the basic con- 
stants of nature are slowly varying functions of the cosmic expansion, as suggested by numerous 
experiments. We have connected the variation of the nuclear and particle masses, fundamental 
scales and particle physics couplings (e.g. the fine structure and the strong coupling "constant") 
to the possible cosmic evolution of the two parameters p\ and G of Einstein's gravity theory, 
i.e. the vacuum energy density (or cosmological "constant") and the gravity constant. The non- 
conservation of matter, associated to a time variation of the parameters in nuclear and particle 
physics, must be compensated by the corresponding evolution of the vacuum energy density and/or 
gravitational coupling G. This resulting picture of the cosmic evolution is compatible with the cos- 
mological principle, but p\ and G evolve with the cosmic time in combination with the fundamental 
"constants" 07]. 

We have represented the possible time evolution of the physical quantities in terms of the 
dimensionless parameter f c fj (proportional to the original v\ If the experiments would detect a 
mass density anomaly in the microphysics world, e.g. through a (red)shift in the value of the 
proton to electron mass ratio, it would lead to a non-vanishing value of vb (which is the baryonic 
part of v e s). This anomaly would be correlated with the corresponding shift of the dimensionless 
quantities A/nip and Gm^ (for the class of Models II and III respectively). A shift in the value of 
these dimensionless quantities would determine v cx f c g, and the corresponding value of v c r could 
be confronted with a possible mass anomaly ub of the nuclear matter. From the difference with 
f e ff we could infer an indirect effect from dark matter, which is controlled by the dimensionless 
index vx- 

We have described the cosmic evolution of the various quantities through the Hubble rate H 
as the basic scale, which can parameterize the running of the masses and couplings as well as the 
vacuum energy and/or Newton's constant G. The running of p\ and G is related to the quantum 
effects of the particles on the effective action of QFT in curved space-time. The vacuum energy 
density is written as a function of H: p\ = p\(H). 

Matter is non-conserved. We have attributed the non-conservation to a cosmic redshift (hence 
a cosmic time variation) of the masses of the nucleons, due to the corresponding change of the 
QCD parameter Aqcd- All atoms would be affected as well. One may expect that the redshift 
should affect the masses of all the fundamental particles (quarks and leptons), including the dark 
matter particles. We have explicitly proposed a connection of the cosmic time evolution of the 
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Aqcd scale and of the elementary particle masses to the corresponding running of p\ and/or G. 

The present bounds, obtained for the time variation of the fundamental constants of nuclear 
matter, point to a rate of change of the nucleon mass and the Aqcd parameter, which is compatible 
with the corresponding bounds on the cosmic evolution of p\ and G. The relevant dimensionless 
parameter, which controls the running of these quantities, must be of order \u\ < 10 -3 or less. The 
current time variation of the vacuum energy can be of order \p\/ pa\ ~ 10~ 3 Hq ~ 10~ 14 yr _1 . This 
can be compared with the current measured rate of change of the Aqcd scale in astrophysical and 
in atomic clock experiments, which provide bounds of the same order of magnitude. However the 
laboratory bounds affect only the nuclear matter contribution to v. The remaining contribution, as 
indicated before, should come from dark matter particles. This approach could eventually provide 
an indirect evidence, that dark matter particles exist. 

The models of the cosmic evolution, discussed in this paper, offer an interesting perspective to 
unify the microphysical and the macrophysical laws of nature. The dark energy is the dynamical 
vacuum energy in interaction with matter. If the dark matter would participate in the cosmic red- 
shift, affecting the baryonic matter, there would be an intimate connection between the evolution 
of the dark matter and of the dark energy. 

The ideas presented here can be tested by different kind of experiments. They could help to 
understand the structure and cosmic behavior of ordinary matter as well as to uncover the mysteries 
of dark matter and dark energy. The small cosmic variation of the physical "constants" may signal 
a connection between the large scale structure of the universe and the quantum phenomena in the 
microcosmos. 
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